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I review recent results in QCD at high energy, emphasizing the role of higher-order computations, power 
corrections, and Monte Carlo simulations in the study of a few discrepancies between data and perturbative 
predictions, and discussing future prospects. 


1. INTRODUCTION 

After more than 25 years of considerable the¬ 
oretical and experimental efforts, it appears that 
QCD is the theory of strong interactions. Ideally, 
in high-energy QCD one needs one single piece 
of information from the experiments: the value 
of a s . Starting from that measured value, ev¬ 
ery observable can be computed from first princi¬ 
ples. In practice this is not feasible, since we don’t 
know how to perform calculations in terms of the 
hadrons that experiments measure in their detec¬ 
tors. Perturbation theory offers a viable way out, 
since it allows to prove, at least formally, the so- 
called factorization theorems. These give explicit 
prescriptions to write physical observables as the 
convolution of short- and long-distance parts, up 
to terms suppressed by the power of some large 
scale. We can imagine this factorization to oc¬ 
cur at an arbitrary scale /x; with a suitable choice 
of n the short-distance pieces, which are entirely 
expressed in terms of quarks and gluons, are per- 
turbatively calculable. The long-distance pieces 
(such as parton densities) cannot be computed 
in perturbation theory, but their dependence on 
/x can. Furthermore, they are universal, which 
means that they don’t depend on short-distance 
physics, but solely on the nature of the hadrons 
involved, which is a key factor for perturbative 
QCD to have predictive power. 

Pending a general solution of QCD, the com¬ 
puting framework based on perturbation theory 
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may be regarded as a hypothesis, which needs to 
be supported, or disproved, by experimental ob¬ 
servations. Countless tests have indeed been suc¬ 
cessful, convincing us of the correctness of this 
approach and of the capability of QCD to de¬ 
scribe strong interactions; in many areas precise 
measurements, rather than tests, are being car¬ 
ried out. This success may give to the non-expert 
the impression that current efforts in theoretical 
QCD are perhaps technically appealing, but not 
compelling physics-wise. To counter this view, it 
is worth reminding that in the past decade the 
studies of several technically-involved problems, 
such as computations to next-to-leading order 
accuracy, resummation techniques, and Monte 
Carlo simulations, have been key factors to the 
outstanding achievements of LEP, SLC, HERA, 
and Tevatron. However, the solutions devised so 
far are not sufficient any longer. For an improve¬ 
ment of the accuracy in the extraction of a s , for 
a deeper understanding of the interplay between 
perturbative and non-perturbative physics, and 
for a realistic modelling of Tevatron Run II and 
LHC physics, new ideas and computations are 
necessary. It must be clear that such investiga¬ 
tions are not only relevant to the study of QCD 
itself, but also to a variety of other issues, from 
SM precision tests to searches of beyond-the-SM 
physics. Besides, a few unsatisfactory results re¬ 
main in QCD, which deserve further studies. 

Needless to say, this review cannot be complete, 
and I’ll have to leave out several interesting re¬ 
sults, such as new NLO calculations, progress in 
small -x physics, diffraction, fully numerical com- 
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putations, and spin physics. Also, I’ll not present 
most of the recent experimental results in high- 
energy QCD, since they can be found in K. Long’s 
writeup [ 0. I’d rather use a few phenomenologi¬ 
cal examples to discuss some theoretical advance¬ 
ments, and open problems. Related topics can be 
found in Z. Bern’s writeup [ 0], I’ll quote pa¬ 
pers submitted to this conference as [S-NNN], S 
and NNN being the session and paper numbers 
respectively. 

2. HEAVY FLAVOURS 

Heavy flavour production is one of the most 
extensively studied topics in QCD. An impres¬ 
sive amount of data is available, for basically all 
kinds of colliding particles. The non-vanishing 
quark mass allows the definition of open-heavy- 
quark cross sections (whereas for light quarks 
one must convolute the short-distance cross sec¬ 
tions with fragmentation functions, in order to 
cancel final-state collinear divergences). On the 
other hand, the presence of the mass makes 
the calculation of the matrix elements more in¬ 
volved. The breakthrough was the computation 
of total QQ hadroproduction rates to NLO ac¬ 
curacy [ |], [|] , readily extended to other produc¬ 
tion processes and more exclusive final states [ 
|, §, 0, § g, 0 [li], [b|. The resummation 
of various classes of large logarithms affecting 
these fixed-order computations, such as thresh¬ 
old, large-p T , and small-x logs, has also been ac¬ 
complished, typically at the next-to-leading log 
(NLL) accuracy. 

In those kinematical regions not affected by 
large logs, the mass of the heavy quark sets the 
hard scale. Furthermore, the impact of effects of 
non-perturbative origin (such as colour drag or 
intrinsic fc T ) is known to be larger the smaller the 
quark mass and CM energy. Thus, top physics is 
expected to be the ideal testing ground for per¬ 
turbative computations. The agreement between 
NLO results [ 0] (dashed lines - the band is the 
spread of the prediction due to scale variation), 
and Tevatron Run I data [H0 - shown in fig. | 
for total tt rates, appears in fact to be satisfac¬ 
tory. The inclusion of soft-gluon effects (solid 
lines), resunnned to NLL accuracy according to 



Figure 1. Total tt rate at the Tevatron. 


the computation of ref. [ [l5|] , is seen to increase 
only marginally the NLO prediction, while siz- 
ably reducing the scale uncertainty. Top produc¬ 
tion appears therefore under perturbative control. 
More stringent tests will be performed in Run II: 
the errors on mass and rate will be smaller, and 
measurements will be performed of more exclu¬ 
sive tt observables and of single-top cross section 
(for which fully differential NLO results are now 
available [ 0 )- 

Bottom quarks are copiously produced at col¬ 
liders, and precise data for single-inclusive distri¬ 
butions have been available for a long time. It 
is well known (see ref. [ 17 for a review) that 
NLO predictions are about a factor of two smaller 
than data at the SpS and at the Tevatron (on 
the other hand, the shape of the p T spectrum of 
the centrally-produced 6 is fairly well described 
by QCD). In a recently-published measurement [ 
0 of the B + p T spectrum, CDF find that the 
average data/theory ratio is 2.9 ± 0.2 ± 0.4. How¬ 
ever, this worrisome result is largely due to an 
improper computation of the NLO cross section. 
Let me remind that the spectrum of a 6-flavoured 
meson B is computed as follows: 

= [ dzdp T D(z; e)-jr-S(p T — zp T ), (1) 
dp T j dp T 


where p T (p T ) is the transverse momentum of b 
(R), dab is the cross section for open-6 produc¬ 
tion, and D{z ; e) is the non-perturbative frag¬ 
mentation function (NPFF), which describes the 
b —> B fragmentation. NPFF is not calcula- 
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ble from first principles, and the free parameter 
it contains (e) is fitted to data after assuming 
a functional form in z (such as Peterson [ 0, 
Kartvelishvili [ |2(| , etc). This fit is typically per¬ 
formed using eq. ( 0 ), identifying the l.h.s. with 
e + e _ data. It follows that the value of e is strictly 
correlated to the short-distance cross section dab 
used in the fitting procedure, and thus is non¬ 
physical When eq. (|J) is used to predict 13-meson 
cross sections, it is therefore mandatory to make 
consistent choices for e and dab■ This has not 
been done in ref. [ |u|: for dab, the NLO result of 
ref. [||] is used, but the value of e adopted (0.006) 
has been derived in the context of a LO, rather 
than NLO, computation. On the other hand, if 
a more appropriate value of e is chosen (~0.002 
[HI), the theoretical prediction increases by a 
mere 20% [ , still rather far from the data. 

There are, however, a couple of observations 
which save the day. First, one has to remark 
that in the upper end of the p T range probed 
by CDF (p T ~ 20 GeV), large-logp T /m effects 
may start to show up. Therefore, FONLL com¬ 
putations [ ^|| should be used rather than NLO 
ones. The FONLL formalism consistently com¬ 
bines (i.e., avoids overcounting) the NLO result 
with the cross section in which p T /m logs are re¬ 
summed to NLL accuracy (such resummed cross 
section is sometimes incorrectly referred to as 
“massless”). Thus, FONLL can describe both the 
small-p T (p T ~ to, where resummed results don’t 
make sense) and the large-p T (p T to, where 
NLO results are not reliable) regimes. The second 
observation concerns again the NPFF: dab / dp T is 
a rather steeply falling function, and one can ap¬ 
proximate it with C/pN in the whole p T range; 
then (from eq. ( 0 )) das/dp T = DjsrC/p^ , where 
Dn = J dzz N ~ 1 D(z ) is the N th Mellin moment 
of the NPFF. This fact has been noticed some 
time ago [ 0 , and DnC/p^ is seen to approx¬ 
imate the exact result fairly well [ Q. Since 
N = 3-5 (at the Tevatron), it follows that, in 
order to have an accurate prediction for the p T 
spectrum in lradroproduction, it is mandatory 
that the first few Mellin moments computed with 
D(z) agreed with those measured. In ref. [ HI , 
it is pointed out that this is not the case, in spite 
of the fact that the prediction for the inclusive b 


cross section in e + e _ collisions, obtained with the 
same D(z), displays an excellent agreement with 
the data. There may seem to be a contradic¬ 
tion in this statement: if the shape is reproduced 
well, why this is not true for Mellin moments? 
The reason is that when fitting D(z) one excludes 
the region of large 2 , since it is affected by Su- 
dakov logs, and by complex non-perturbative ef¬ 
fects which are unlikely to be described by the 
NPFF. On the other hand, the large -2 region is 
important for the computation of Djv (because 
of the factor 2 JV_1 in the integrand). Therefore, 
for the purpose of predicting 33-meson spectra at 
colliders, ref. [ |22| advocates the procedure of fit¬ 
ting the NPFF directly in the iV-space. A fit to 
the second moment (denoted as N = 2 fit hence¬ 
forth) is found to fit well all the D^s for N up 
to 10; and although Kartvelishvili’s form is used, 
Peterson’s gives comparable results. Using the 



Figure 2. B + data vs theory 



FONLL computation, and a N = 2 fit for the 
NPFF, the average data/theory ratio reduces to 
1.7 ± 0.5 ± 0.5 [ H^]. Taking the scale uncertainty 
into account, B + data appear to be compatible 
with QCD predictions (see fig. |[). 

If one wants to avoid the pitfalls of NPFF’s, an 
alternative possibility consists in considering b- 
jets rather than B mesons, since in this case the 
NPFF simply doesn’t enter the cross section. The 
comparison between NLO predictions for &-jets [ 
p5| and DO measurements [ |2(| is indeed satisfac¬ 
tory: data are consistent with theory in the range 
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25 < E^~-’ et < 100 GeV. Overall, one can con¬ 
clude that b data at the Tevatron are reasonably 
described by NLO QCD. It is worth mentioning 
that some existing results, presented in terms of 
6 -quark cross sections, are likely affected by the 
findings of ref. [ 0 ] , and need to be reconsidered. 
Among the various mechanisms which can fur¬ 
ther increase the theoretical predictions, small-a: 
[ |28| and threshold resummations [® will 

probably play a secondary role wrt NNLO contri¬ 
butions, which are expected to be large given the 
size of the K-factor at the NLO. 

I now turn to the case of charm production. 
A thorough discussion on this topic is beyond the 
scope of this review, and I’ll only give the briefest 
of the summaries (which will not do any justice to 
the field). LEP data for total rates are in agree¬ 
ment with NLO QCD predictions [ [?J ; the shapes 
of single-inclusive D* spectra in 77 collisions are 
as predicted by NLO QCD [ |||], whereas nor¬ 
malization is off by a factor 1.5-2, but still con¬ 
sistent with QCD when theoretical uncertainties 
are taken into account. The vast majority of fixed 
target hadro- and photoproduction data are well 
described by NLO computations, but only if pre¬ 
dictions for single-inclusive distributions and cor¬ 
relations are supplemented by some parametriza- 
tions of non-perturbative phenomena (such as in¬ 
trinsic k T ). At HERA, DIS data are in agreement 
with NLO QCD results [ ^2[. In photoproduction, 
some concerns have arisen in the past because of 
the discrepancy between ZEUS and HI measure¬ 
ments in the comparison with theory: HI [ |30| 
appears to be in perfect agreement with QCD, 


whereas ZEUS [ 31 is at places (for intermedi¬ 
ate p T ’s and large rj's) incompatible with NLO 
predictions. ZEUS have submitted to this confer¬ 
ence [5-786] data with unprecedented coverage at 
large p T . The comparison to FONLL predictions [ 
0 , 0| , shown in fig. |j, appears to be satisfactory, 
although data are marginally harder than theory. 
The agreement improves if a A” = 2 fit for the 
NPFF is adopted (this result is still preliminary). 


Let me finally mention the increasing amount 
of measurements for b rates from fixed-target [ 
|J, HERA [ 0 0 [5-783,5-784,5-785,5-1013,5- 
1014], and LEP [ |37|, [5-366,5-475] experiments. 



Figure 3. ZEUS D* data in 7 p collisions vs 
FONLL predictions 


A summary of the situation, in the form of ratios 
data/NLO QCD, is presented in fig. ||. While the 
fixed-target measurements are in overall agree¬ 
ment with QCD, HERA and LEP measurements 
are largely incompatible with theory. I find this 
hard to reconcile with the results presented so 
far, and the size of the discrepancy also makes 
it problematic to find an explanation in terms of 
beyond-the-SM physics, let alone higher orders in 
QCD. It is necessary to note that in many cases 



Figure 4. Ratios data/theory for b rates. 


the experimental results are extrapolated to the 
full phase space from a rather narrow visible re¬ 
gion. It is encouraging to note that in the cases 
of the recent ZEUS measurements in DIS [5-783] 
and photoproduction [5-785] (full boxes in fig. [|), 
results for single-inclusive distributions are given 
too, which are seen to be fully compatible with 
the corresponding NLO predictions. 
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3. POWER CORRECTIONS 

The necessity of understanding long-distance 
effects in final-state measurements is not peculiar 
to b hadroproduction. The hadron-parton dual¬ 
ity assumption states that there is a class of ob¬ 
servables (such as jet variables or event shapes 
in e + e~ collisions) whose description in terms 
of quarks and gluons is expected to reproduce 
the data, up to terms suppressed by some in¬ 
verse power of the hard scale Q of the process 
(power corrections). These terms are usually es¬ 
timated by comparing the parton- and hadron- 
level predictions of Monte Carlo (MC) generators. 
This procedure is not really satisfactory, since 
MC parameters are tuned to data (which creates 
a bias on the “predictions” for power-suppressed 
effects), and since the definition of parton- and 
hadron-level is far from being straightforward. 

It is remarkable that we can get insight on non- 
perturbative physics from perturbative considera¬ 
tions. The perturbative series, being asymptotic, 
can be summed to all orders only after defining a 
summation procedure (in a rather arbitrary man¬ 
ner); one assumes that this technical manipula¬ 
tion mimics the role played in Nature by non- 
perturbative effects, which are necessary for QCD 
to be self-consistent. The summation procedure 
must eliminate the divergence of the perturbative 
series, but some finite quantities are left uncon¬ 
strained. Thus, the idea is to use the ambigu¬ 
ities of the summation procedure to study non- 
perturbative effects. Although the regularization 
of the divergence can be technically very compli¬ 
cated, it can always be seen as a prescription to 
deal with the Landau pole of a s . The idea of 
ref. [ (DMW from now on) is to bypass 

such a prescription by defining a s in the infrared 
(IR) region, assuming its universality; thus, a s 
should effectively measure confinement effects in 
inclusive quantities (in order to use such an a s 
in actual computations, it is also necessary to 
assume that the concepts of quarks and gluons 
still make sense in the IR). After giving the gluon 
a fake (“trigger”) mass p,, the (fully inclusive) 
observable under study is computed in pertur¬ 
bation theory; the small-/u behaviour determines 
the power p of the leading power-suppressed term 


A p /Q p . The coefficient A p cannot be computed, 
but can be expressed as an integral of k p ~ 1 a s (k) 
over 0 < k < Pi, with /q ~ O(GeV). Since all 
the (non computable) power-correction effects are 
contained in this coefficient, one effectively gets a 
factorization formula. A class of observables of 
great physical relevance is that of event shapes 
for which p = 1; their mean values can be com¬ 
puted with the DMW approach. The property of 
a s universality is formally even more far-reaching; 
for example, it implies that the non-perturbative 
effects it describes exponentiate for observables 
4f| . This offers the possibility 
of studying not only mean values, but distribu¬ 
tions. The results of DMW for mean values can 
then be recovered by an expansion of the (Su- 
dakov) exponent and subsequent average. More 
precisely, if T denotes the observable, the Su- 
dakov is expanded in the region q/Q < T C 1, 
and only the first non-trivial term in the expan¬ 
sion is kept. 

In this context, factorization derives from the 
hypothesis of a s IR universality, which is an ef¬ 
fective description of long-distance effects, and as 
such must be insensitive to the details of parton 
dynamics. For example, these details are irrel¬ 
evant when one sums inclusively over the decay 
products of any parton branchings. It has been 
observed [ ||| that inclusiveness is lost to a certain 
extent when recoil effects (i.e., higher orders) are 
considered, and thus factorization breaks down in 
this case. One may still insist that factorization 
holds, and devise a procedure to systematically 
account for those effects which would break it in 
a naive treatment [ |d|, |b], Q. The results for 
mean values and distributions of event shapes can 
be written as follows: 


that do so [ Eo, 41, 


(fT) — (T)pert + Ct'P , 

(T - ctV), 


dT [ 1 dT 


pert 


P = i£L M t ‘ 1 


Ql 


- a s (Q) + 0(al(Q)) 


/■Mi 

p, I a 0 {pi)= / dka s (k). 

J o 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


Here, cr is a computable coefficient, “pert” 
means perturbatively-computed, and M. includes 
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the two-loop results of refs. [ |d|, |I|, in 
principle, M. can depend upon T, but (acciden¬ 
tally) it does not. The a s (Q) and a^(Q) terms 
in eq. (|]) are there to avoid double counting with 
the perturbatively-computed part, which shows 
again that long- and short-distance effects are cor¬ 
related. Since no small parameter is involved in 
the two-loop computation of At, one may won¬ 
der whether factorization could be spoiled be¬ 
yond two loops. Although it is argued that this is 
not the case [ Q , a comparison with the data is 
mandatory. Also notice that the formulae above 
need to be modified when more complicated kine- 
matical effects have to be described, as in the case 
of broadenings [ !!• Eqs. (|j) ([]]) are used to 
fit the data in terms of a s (Q) and ao(ni). Up¬ 
dated analyses relevant to e + e _ collisions have 
been presented to this conference [5-228,5-229,5- 
389], [ Q - see also [ fiof . Results obtained 
from mean values are satisfactory, with compara¬ 
ble op's obtained from different observables, and 
a s (Mz) values fairly consistent with the world 
average. The situation worsens in the case of dis¬ 
tributions: op universality holds at ~25% level 
(1-2 er), and a s (Mz) values are systematically 
lower than the world average, especially for ob¬ 
servables such as wide broadening Bw and heavy 
jet mass Mh■ This fact is disturbing since the 
same data for distributions, with hadronization 
effects described by MC’s, return a s (Mz) values 
in much better agreement with the world aver¬ 
age. These findings are summarized in fig. |^, 
where the results obtained with the DMW ap¬ 
proach are shown as boxes and crosses for mean 
values and distributions respectively; the world 
average a s (Mz) = 0.1184 ± 0.0031 [ ||(J is also 
shown. Data have been taken in the PETRA, 
LEP and LEP2 energy range, and analyses have 
been presented by Aleph [5-296], Delphi [5-228,5- 
229], Jade [5-389] [ |§, L3 [5-495], and Opal [5- 
368]. One of the two results presented in ref. [ 
|4S|| , reported as the uppermost cross in fig. i is 
obtained by excluding Bw from the fit. Recently, 
the NLL resummation of many DIS event shapes 
has been achieved ([ ^jj, and references therein). 
Power-correction effects can then be studied simi¬ 
larly to what done for e + e - collisions, and the re¬ 
sults are intriguing. Fig. (taken from ref. [ HI) 



°<s(M z ) 

Figure 5. Results for a s (Mz) from event shapes 
in e + e^ collisions. 


presents the a s (Mz) and op values obtained from 
fitting event shape distributions in DIS, and event 
shape means in e + e _ . The DIS and e + e~ results 
are largely consistent, which implies, in view of 
what shown in fig. |^, that event shape distribu¬ 
tions in e + e _ and DIS prefer different values for 

ot s {M z )- 



Figure 6. Results for a s (Mz) and op from DIS 
and e + e _ event shapes. From [ |52f. 


Although not compelling from the statistical 
point of view, these results for event shape dis- 
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tributions in the DMW approach may hint to 
the necessity of a more complete description 
of hadronization effects. As mentioned before, 
eq. (|ij) results from keeping the first non-trivial 
term in a Taylor expansion. If no expansion 
is made, from rather general factorization argu¬ 
ments in the two-jet limit T —> 0 one gets the 
following formula [ |Ai|| 


— (T) = 
dT y ’ 


rTQ 


d eMe) * 


(T -e/Q), ( 6 ) 


pert 


where fr{z) is known as shape function. DMW 
formulae are recovered with fr (e) = 5(e—QctV ); 
in the general case, the first Mellin moment of 
fr has the same meaning as ao of DMW. With 
eq. ®) it is not necessary to assume that TQ>/Ji 
(all terms 1 /(TQ) n are now expected to be in¬ 
cluded), and therefore the fit ranges can be ex¬ 
tended. Similarly to ao, fr cannot be computed 
from first principles; thus, in order not to lose 
predictive power, a functional form depending on 
a small number of parameters must be assumed [ 
pI~T[ . [55| , p| , keeping in mind that QCD dynamics 
and Lorentz invariance considerations [[57j can be 
used to severely constrain the form of a more gen¬ 
eral shape function, from which fr is derived, in¬ 
dependently of phenomenological arguments. In 
a different approach (DGE [ 55, 58) ), whose final 
result has the same form as eq. ®), it is sug¬ 
gested to combine Sudakov and renormalon re- 
summations in a single formalism. A renormalon 
ambiguity appears in the exponent, and the pre¬ 
scription necessary to resolve it can be naturally 
formulated in terms of a shape function, auto¬ 
matically constraining its functional form. Al¬ 
though the fr which one gets in DGE is con¬ 
sistent with the one obtained in refs. 0 0 - 
it has to be stressed that the perturbative result 
da / dT\ per t in eq. (p) is different in the two for¬ 
malisms, since DGE includes a class of sublead¬ 
ing logs of renormalon origin. In e + e _ physics, 
eq. (^j) gives satisfactory results: a good fit to 
the second moments of 1 -thrust, Mjj, and C pa¬ 
rameter is obtained in ref. [ p| , and the fits for 
thrust and Mh of ref. [ p| are in better agree¬ 
ment in the (a s ,ao) plane wrt those obtained 
with DMW. On the other hand, according to ref. 
[ |1, a s (M z ) = 0.1086 ± 0.0004(ea:p) (the the¬ 


ory error is estimated to be around 5%). There¬ 
fore, it seems that a more refined treatment of 
non-perturbative effects, which is helpful in other 
respects, is not what one needs in order to get 
larger a s values. 

In a couple of interesting analyses, Delphi 
adopted rather unconventional methods to study 
event shapes. In ref. [ ||ij, event shape dis¬ 
tributions were compared to fixed-order 0(cq) 
(NLO) results (i.e., resummation has not been in¬ 
cluded), using the renormalization scale as a free 
parameter in the fit, and correcting for hadroniza¬ 
tion effects with MC’s. Although I’m not aware 
of any theoretical consideration which justifies 
such a procedure (called in ref. [ 0 “experi¬ 
mental optimization of the scale”), the a s val¬ 
ues obtained from different observables display a 
remarkable consistency, and they are also consis¬ 
tent with those obtained by using NLL-resummed 
predictions. In [5-228] (see also ref. [ @) event 
shape means were shown to give mutually consis¬ 
tent a s values in the context of a renormalization 
group approach (RGI [ 0), without needing any 
hadronization corrections. The final a s results of 
refs. [ |5^] and [5-228] are in excellent agreement 
with the world average. I interpret these findings 
as the indication that (at least in a given CM en¬ 
ergy range) the uncertainties affecting theoretical 
predictions at the NLO are larger than or of the 
same order as the power-suppressed effects that 
one aims to study. It thus appears that the com¬ 
putation of event shapes at the NNLO is neces¬ 
sary for a deeper understanding of this matter. 

In summary, in the past few years a solid 
progress has been achieved in the understand¬ 
ing of power-suppressed effects in e + e _ collisions 
and in DIS. Although models such as DMW 
can successfully describe many features of the 
data, some aspects deserve further studies. In 
some cases, improvements are obtained within ap¬ 
proaches which refine the treatment of the non- 
perturbative part, using a shape function, but the 
computation of the next order in perturbation 
theory will likely be necessary in order to obtain a 
more consistent overall picture. It is worth recall¬ 
ing that the study of hadron-mass effects has been 
found [ ^ 2 | to induce further power-suppressed 
terms, some of which can be eliminated by adopt- 





ing a suitable definition for the observables (E- 
scherne). Furthermore, more stringent tests of 
the models for power-suppressed effects should 
be performed using observables with more com¬ 
plicated kinematic structure and/or gluons at the 
LO (see refs. [0, 0 and references therein). Fi¬ 
nally, it appears to be mandatory to extend the 
studies of such models to the case of hadronic col¬ 
lisions (for jet observables in particular), where 
little work has been done so far. 


4. NNLO COMPUTATIONS 


Bottom production at the Tevatron and event 
shapes in e + e - collisions are a couple of exam¬ 
ples which provide physical motivations to in¬ 
crease the precision of the perturbative compu¬ 
tations. If tree n-point functions contribute to a 
given reaction at the LO, the N fe LO result (i.e., of 
relative order a* wrt to the LO) will get contri¬ 
butions from the l- loop, (n + p)-point functions, 
with l + p < k. There are basically three ma¬ 
jor steps to make in order to get physical predic¬ 
tions: i) explicit computation of all the /-loop, 
(n + p)-point functions; ii ) cancellation of soft 
and collinear divergences (which I’ll denote - im¬ 
properly - as IR cancellation henceforth); Hi) nu¬ 
merical integration of the finite result obtained 
from i) and ii), with MC techniques to allow 
more flexibility. One should also mention that 
UV renormalization has in general to be carried 
out; however, this is basically textbook matter by 
now, and thus I’ll not deal with it in the follow¬ 
ing. For NLO computations (k = 1), steps i)- 
iii) appear to be understood. One-loop integrals 
have been computed up to five external legs [ 0 ) ; 
the case of n > 6 cannot probably be handled 
with Feynman-diagram techniques only, and still 
awaits for a general solution (see [ ||). Subtrac¬ 
tion [ 0) and slicing [ 67 methods, to achieve 
IR cancellation with semi-analytical techniques, 
have been available for a long time. In their mod¬ 
ern versions 0 - 00 - 0 ^ 0 0 - 0 they are 
formulated in an universal (i.e., process- and n- 
independent) way, which simplifies step Hi) con¬ 
siderably, and allows the computation of any IR 
safe observable (no matter how exclusive). No 
modification is needed in order to incorporate 


new one-loop results. Attempts to achieve IR 
cancellation through full numerical computations 
[ 0 - 0 , 0 are still in a preliminary stage, re¬ 
producing known results for three-jet production 
in e + e^ collisions. 

Only a handful of production processes have 
been computed to NNLO accuracy and beyond: 
results are available for DIS coefficient functions 
[ 0 0 and for the Drell-Yan K-factor [0, 81] 
at 0(a g), and for the rate e + e _ —> hadrons [ 82] 
at O(af) (I’ll later deal with inclusive Higgs pro¬ 
duction in some details). All these computations 
are inclusive enough to allow a complete analyti¬ 
cal integration over the phase-space of final-state 
partons. Such an integration is not possible in 
general, and an NNLO process-independent and 
exclusive formulation of IR cancellation will likely 
proceed through semi-analytical techniques, sim¬ 
ilar to those adopted at the NLO. Therefore, one 
needs to know the IR-divergent pieces of all the 
quantities contributing to the cross section at the 
NNLO. These are (I still assume that the LO gets 
contribution from the tree n-point functions): 
a) tree-level (n + 2 )-parton amplitudes squared 
|T„ + 2 | ; b) the interference between tree-level 
[n + l)-parton amplitudes and one-loop [n + 1 )- 
parton amplitudes 3?(T* + 1 l| 1 1 | 1 ); c) the interfer¬ 
ence between tree-level n-parton amplitudes and 
two-loop n-parton amplitudes 3ff (T*L y n ’)-, d) one- 
loop n-parton amplitudes squared |L^| 2 . The 
IR divergences of |T n+ 2| 2 result from having two 
soft partons, or three collinear partons, or one 
soft parton plus two other collinear partons, or 
two pairs of collinear partons. Only the latter 
configuration is trivial, in the sense that the cor¬ 
responding singular behaviour of |T n+ 2| 2 can be 
obtained from known NLO results; the other lim¬ 
its have been studied in refs. [ 0 , 0 , 0 , 0 0 ] , 
and can be generally cast in the form of a re¬ 
duced n-resolved-parton matrix element, times a 
suitable kernel. The problem of combining these 
singular pieces into local IR counterterms, and of 
integrating these counterterms over the appropri¬ 
ate region of the phase space, is still unsolved. 
One also needs to know the IR divergences of 
when one parton is soft or two partons are 
collinear: this is a new feature of NNLO cornpu- 
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tations, since at the NLO all partons in a virtual 
contribution are resolved. These limits are also 
known [ [s8[ [dj| [t]. [Of], |9^, ||| [H]| . Finally, the 
general form of the residues of the poles 1 /e 4 , 
1 /e 3 , and 1 /e 2 appearing in Ln' 1 has been given 
in ref. [ |95| , without computing any two-loop in¬ 
tegrals (see also [ @). 

The pole terms found in ref. [ |9^] must pre¬ 
cisely match those resulting from the explicit 
computations of two-loop integrals. A lot of 
progress has been made in the past couple of years 
in such computations, and now all the two-loop 
2 —> 2 and 1* —> 3 amplitudes are available. First, 
all the (very many) tensor integrals are reduced to 
a much smaller number of master scalar integrals, 
thanks to integration-by-part identities [ |||, 
(previously used for two-point functions), and 
Lorentz-invariance identities [ |9^|; integration- 
by-part ident ities for n-leg, /-loop integrals were 
also shown [ | 1 Q[| to be equivalent to those for 
(n — to)- leg, (l + m)-loop integrals. The problem 
of actually computing the master integr als i s of a 
different nature. The breakthrough [ 101, 102 


was the use of a Mellin-Barnes representation 
for the propagators in the computation of pla¬ 
nar and non-planar massless double box integrals 
(expanded in the dimensional-regularization pa¬ 
rameter e). Negative space-dimension techniques 
can also be applied to simpler topologies [ 
1051. The computation of master integrals 


103 


104 


is mapped onto the problem of solving differen¬ 
tial equations in the approach of refs. [ 106, 107|. 


This approach has been adopted to compute all of 
the double box master integrals with one off-shell 
leg 


|104 109 1 , not all of which had been com¬ 
puted with Mellin-Barnes techniques [ 110 , 111 |. 
The master integrals so far computed, together 
with the reduction-to-master-integral techniques, 
would allow the computation of e + e“ —> 3 jets, 
and of two-jet production in hadronic collisions 
( and a few other processes: see for example ref. 
[ 1121 for a discussion - unfortunately, these pro¬ 
cesses don’t include heavy flavour hadroproduc- 
tion) if one knew how to achieve IR cancellation 
for a generic observable at the NNLO. 

In hadronic physics, the computation of NNLO 
short-distance cross sections is not sufficient to 


get NNLO-accurate predictions, since NNLO- 
evolved PDF’s are also necessary. NNLO-evolved 
PDF’s require the computation of Altarelli-Parisi 
splitting functions to three loops. It turns out to 
be convenient to perform such a computation in 
Mellin space; the results for the first few Mellin 
moments [ 113, 114, |l!5|| (together with con¬ 


straints on the small-a; behaviour) have been used 
[ 116, 117, 118{ to obtain approximate expressions 
for the splitting functions in the x space. Very 
recently, the complete three-loop computation of 
the n F part of the non-singlet structure function 


in DIS has become available [ 119 , from which the 
corresponding coefficient functions (relevant to 
N 3 LO computations) and splitting function can 
be extracted; the latter has been cross-checked 
against the approximate results mentioned above, 
and full agreement has been found. Although the 
complete expressions fo r th e splitting functions 
will not appear soon [ |ll9| , this fairly impres¬ 
sive result gives confidence on t he a ccuracy of the 
approximate solution of ref. [ |l! 8 | . One has to 
keep in mind that, in order for any PDF set (such 
as NNLO-MRST [ 1201) to actually be of NNLO 


accuracy, not only the three-loop splitting func¬ 
tions are needed, but also all the short-distance 
cross sections used in the fits must be computed 
to NNLO. At present, this is the case for the DIS 
coefficient functions only (which implies that the 
approximation is generally good, the bulk of the 
data being from DIS); for example, in the case of 
Drell-Yan the x F distribution needed for the fits 
has only NLO accuracy (the Drell-Yan NNLO Re¬ 
factor is used for normalization purposes [ 120 |). 


This is another hint of the necessity of solving 
the problem of IR cancellation at the NNLO in a 
general way. 

In view of enormous amount of work done and 
yet to be done, one may ask whether the final 
outcome will justify such an effort. The answer 
is certainly positive, but one must keep in mind 
that NNLO computations will not automatically 
mean precision physics. As shown in the case 
of e + e _ event shapes, short- and large-distance 
effects are always correlated, and any advance¬ 
ment in perturbation theory should be comple¬ 
mented by a deeper understanding of this correla¬ 
tion. Furthermore, no precision study in hadronic 
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collisions can be made without an accurate assess¬ 
ment of uncertainties due to PDFs. This matter 
is now receiving considerable attention: see for 
example ref. [ 112] for a review. NNLO compu¬ 
tations will certainly play a major role in those 
cases in which NLO results still give an unsat¬ 
isfactory description of data: examples are the 
processes with large K-factors (such as b produc¬ 
tion), or the observables which require a better 
description in terms of kinematics (such as jet 
profiles). 

A large-K-factor process is the direct produc¬ 
tion of SM Higgs at hadron colliders. At the NLO, 
the exact computation [ 1211 for the gg channel is 
found to be i n ex c ellen t agreement with approxi¬ 
mate results [ 122, 123 based on keeping the lead¬ 
ing term of an expansion in m,,/m;r op (the agree¬ 
ment further improves if the full dependence on 
m\/m% op is kept in the LO term). Therefore, one 
assumes the same to hold at the NNLO, and com¬ 
putes the NNLO contribution in the mt op —> oo 
limit. This is feasible since the effec t ive La- 
grangian ggH has been obtained [ 124 , 125 ] to 
O(af) (actually, one order larger than nece ssary 
here). Exploiting the technique of ref. [ |10(| . the 
two-loop virtual correction to gg ^ H has been 
obtained in ref. [ 1261. The missing contributions 
to the physical cross section (including qg and 
qq channels) have been presented in ref. [ |127| , 
where an expansion for x H ~ 1 {x n = to|/s) has 
been used to compute the double-real contribu¬ 
tion. Terms up to (1 — x H ) 16 have been included. 
The x H — 1 expansion is expected to work well, 
since the gluon density is rather peaked towards 
small x’s , and thus the CM energy available at 
the partonic level is nev er too far from threshold. 
As argued in ref. [ 12S], collinear radiation also 
gives a sizable effect. Me thods sim ilar to [ 125 ] 
have been used in refs. [ 12E, 1301, with the re¬ 


sult of ref. [ |l2q| , to give early estimates of the 
complete NNLO rate. These estima tes a re seen 
to agree well with the result of ref. [ 127|, which 
confirms the soft-collinear dominance in inclusive 
Higgs p rodu ction at colliders. Finally, the result 
of ref. [ 127| has been found to agree to 1% level or 
better with the computation of ref. [ 1311, where 
the double-real contribution is also evaluated ex¬ 
actly (it is interesting to notice that in ref. [131] 


the phase-space integrals relevant to real-emission 
terms are computed with techniques used so far 
only for loop diagrams - it remains to be seen 
whether this method can be generalized to more 
exclusive observables). As shown in fig. the 



Figure 7. SM Higgs total rate at the LHC. 


inclusion of NNLO corrections seems to suggest 
that effects beyond this order are negligible. The 
scale dependence is red uced wrt the one observed 
at the NLO (see ref. [ 131]). The dominance of 
the region x H ~ 1 implies the potential relevance 
of soft-gluon resummation. Preliminary results [ 
|132| , 112] indicate that NNLL resummation en¬ 
hances the NNLO rate by 5-6% (12-15%) at the 
LHC (Tevatron), for 100 < m H < 200 GeV. 

The result for the fully-inclusive rate also serves 
to compute a slightly less inclusive observable, 
namely the rate for Higgs+jets, with the p T of 
any jets imposed to be smaller than a fixed quan¬ 
tity p^ >eto ' 1 (jet veto). Such an observable, which 
should help in reducing the background due to the 
decay ch anne l H —> W*W*, has been computed 
in ref. [ 133] by subtracting the anti-vetoed jet 
cross section p T > p^' eto ^ (obtained in ref. [ 134 ]) 
from the inclusive NNLO result discussed so far. 
The study of more exclusive observables, which 
implies the understanding of IR cancellation at 
NNLO, will certainly prove useful in the future. 
In this case, the dominance of the region x H ~ 1 
might not be as strong as in the case of inclu¬ 
sive rates. Although it is unlikely that the region 
x H ~ 0 will play any role in phenomenological 
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studies, it is worth recalling that in this region 
the approximation m t0 p —»■ oo is not expected to 
work well: in the full theory the dominant contri¬ 
bution for x n —> 0 is single-logarithmic, whereas 
double logs are also found in the large-TOt op the¬ 
ory. The latter terms have been identified explic¬ 
itly in ref. [ 13E] with fc T -factorization arguments; 
at the NNLO, they are seen to coincide with those 
resulting from the explicit computation of ref. [ 
|131|, thus providing a cross check impossible to 


achieve in the comparison of ref [131] with ref. [ 


5. MONTE CARLO SIMULATIONS 

Monte Carlo (MC) programs are essential tools 
in experimental physics, giving fully-fledged de¬ 
scriptions of hadronic final states which cannot be 
obtained in fixed-order computations. Schemat¬ 
ically, an MC works as follows: for a given pro¬ 
cess, which at the LO receives contribution from 
2 —> no reactions, (2 + no)-particle configurations 
are generated, according to exact tree-level ma¬ 
trix element (ME) computations. The quarks and 
gluons (partons henceforth) among these primary 
particles are then allowed to emit more quarks 
and gluons, which are obtained from a parton 
shower or dipole cascade approximation to QCD 
dynamics. This implies that MC’s cannot sim¬ 
ulate the emission of final-state hard (i.e., with 
large relative transverse momenta; thus, hard is 
synonymous of resolved here) partons other than 
the primary ones obtained from ME computa¬ 
tions. Furthermore, total rates are accurate to 
LO. 

Although these problems are always present in 
MC simulations, they become acute when CM 
energies grow large, since in this case channels 
with large numbers of well-separated jets are phe¬ 
nomenologically very important and, correspond¬ 
ingly, total rates need to be computed to an ac¬ 
curacy better than LO. Two strategies can be de¬ 
vised in order to improve MC’s. The first aims at 
having n E extra hard partons in the final state; 
thus, in the example given above, the number of 
final-state hard particles would increase from no 
to no + n E . This approach is usually referred to 
as matrix dement corrections , since the MC must 


use the (2+no+n E )-particle ME’s to generate the 
correct hard kinematics; more details are given 
in sect. 0 The second strategy also aims at 
simulating the production of no + n E hard par¬ 
ticles, but improves the computation of rates as 
well, to N raE LO accuracy. A discussion is given in 


sect. 5.2. 


5.1. Matrix element corrections 

There are basically two major problems in the 
implementation of ME corrections. The first 
problem is that of achieving a fast computation 
of the ME’s themselves for the largest possible 
no + n E , and an efficient phase-space generation. 
The second problem stems from the fact that 
multi-parton ME’s are IR divergent. Clearly, in 
hard-particle configurations IR divergences don’t 
appear; however, the definition of what hard 
means is, to a large extent, arbitrary. In prac¬ 
tice, hardness is achieved by imposing some cuts 
on suitable partonic variables, such as p T ’s and 
(77, (p)-distances cLR in hadronic collisions. I col¬ 
lectively denote these cuts by S sep . One assumes 
that n hard partons will result (after the shower) 
into n jets; but, with a probability depending on 
S sep , a given n-jet event could also result from 
n+m hard partons. This means that, when gener¬ 
ating events at a fixed no + n E number of primary 
particles, physical observables in general depend 
upon 8 S ep ; I refer to this as the <5 sep -bias problem. 
Any solution to the 5 sep -bias problem implies a 
procedure to combine consistently the treatment 
of ME’s with different no + n E s. Here, the diffi¬ 
culty is that of avoiding double counting, that is, 
the generation of the same kinematical configura¬ 
tion more often than prescribed by QCD. 

The vast majority of recent approaches to ME 
corrections address only the first of the two prob¬ 
lems mentioned above. A considerable amount of 
work has been devoted to the coding of hadronic 
processes with vector bosons/Higgs plus heavy 
quarks in the final state, which cannot be found 
(regardless of the number of extra partons) in 
standard MC’s. The complexity of hard-process 
generation for large no + n E suggests to imple¬ 
ment it in a package (which I call ME generator) 
distinct from the shower MC. The ME genera¬ 
tor stores a set of hard configurations in a file 
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(event file); the event file is eventually read by 
the MC, which uses the hard configurations as 
initial conditions for the showers. The advantage 
of this procedure is that it is completely modu¬ 
lar: one given event file can be read by differ¬ 
ent MC’s, and conversely one MC can read event 
files produced by different ME generators. It is 
clearly convenient to reach an agreement on the 
format of such event records: this is now available 


136| (Les Houches accord #1). Ready-to-use 


ME generators (with different numbers of hard 
processes implemented) are AcerMC [ |137| , ALP- 
GEN [ |138| , |l 39| . [14 C| , and MadGraph/MadEvent 
[ |145| , [146| . Related work, at present set up to 
function only wit h Pythia , has been pre sent e d by 
the CompHEP [ 141 , 142 and Grace [ L43, 144] 
groups. All of these ME generators use Feynman- 
diagram techniques in the computation of ME’s, 
except ALPG EN, which uses the iterative algo¬ 
rithm Alpha [ 147 (see also [ 1481). 

When using an ME generator, the cuts S sep 
must be looser than those used to define the ob¬ 
servables. For example, the p T -cut imposed at 
the parton level must be smaller than the min¬ 
imum p T of any jets. On the other hand, the 
cuts should not be too loose: the looser the cuts, 
the larger the probability of getting a n-jet event 
starting from n + m hard partons. Thus, S sep 
must be chosen in a range which is somewhat de¬ 
pendent upon the observables that one wants to 
study. This implies that, strictly speaking, the 
combination of an ME generator with a shower 
MC is not an event generator, since the event 
record depends upon the observables. This hap¬ 
pens precisely because such a combination is af¬ 
fected by the J sep -bias problem: it is therefore 
necessary to assess its impact on physical observ¬ 
ables. An example is given in fig. ||, obtained 
using ALPGEN+Herwig (any other ME genera¬ 
tor and shower MC would give equivalent results 
for the same observable). The plot presents jet 


rates, integrated over E? et ^ > E T o, for the hard¬ 
est jet in W+3-jet events at the Tevatron, versus 
the parton separation dR par t imposed at the level 
of ME generation. Jets are reconstructed with 
the cone algorithm, with R = 0.7. Rates are nor¬ 
malized to the result obtained with dR par t = 0.7. 


The conclusion is that, in the “reasonable” range 
0.3 < dR par t < 0.7, the physical prediction has a 
0(20%) dependence on dR par t■ The existence of 
a 5 sep dependence should always be kept in mind 
when using an ME generator + shower MC com¬ 
bination, because it affects the precision of the 
prediction. On the other hand, this is a rather 
modest price to pay: for multi-jet observables, 
ordinary MC’s can underestimate the cross sec¬ 
tion by orders of magnitude, and the use of ME 
generators is mandatory. 



Figure 8. Dependence of (normalized) jet rates 
on interparton separation. 


As mentioned before, the 5 sep bias can be 
avoided by suitably combining the generation of 
ME’s with differ ent np + n E . Early pro posal for 
ME corrections [ 149, 150 , 151, 152, 1531 achieved 
this in the case n E = 0,1. The solution for ar¬ 
bitrary n E appears to be more complicated; it 
has been fully i mplement ed for shower MC’s in 
e + e _ collisions [ 154, 155] in the case of jet pro¬ 
duction; alon g sim ilar lines, proposal for colour 
dipole MC’s|_156| and shower MC’s in hadronic 
collisions [157] have also been made. The idea 
of ref. [ 1541 is the following, a) Integrate all 
the 7* —> 2 + n E ME’s by imposing yij > y 1N1 
for any pairs of partons i.j. with y im a fixed pa¬ 
rameter and ytj = 2min(E 4 2 , Ej)(l — cos%)/Q 2 
the interparton distance defined according to the 
b) Choose statistically an 


/c T -algorithm [ 158 


n E , using the rates computed in a), c) Generate 
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a (2 + n E )-parton configuration using the exact 
7 * —> 2 + n E ME, and reweight it with a suit¬ 
able combination of Sudakov form factors (corre¬ 
sponding to the probability of no other branch¬ 
ings). d) Use the configuration generated in c) 
as initial condition for a vetoed shower. A vetoed 
shower proceeds as the usual one, except that it 
forbids all branchings i —> jk with yjk > 2 /ini 
without stopping the scale evolution. In ref [ 1541, 
?/i NI plays the role of S sep . Although the selec¬ 
tion of an n E value has a leading-log dependence 
on j/iNi, it can be proved that this dependence 
is cancelled up to next-to-next-to-leading logs in 
physical observables [ 154]. This is illustrated in 
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Figure 9. Dependence of mean D parameter on 

2/i = 2 /ini- 


Mean D parameter 



fig. |9[ wh ere the mean value of the D parame¬ 
ter [ |l59l |scl (full circles) is plotted versus the 
value of 2/ini , and compared to the measurement 
of ref. [ |l6(| (band). Figure || clearly documents 
that the consistent combination of ME’s with dif¬ 
ferent no + n E largely reduces the dependence of 
observables on S sep . This, however, comes at the 
price of modifying the shower algorithm. Further¬ 
more, the procedure is more computing-intensive, 
since all n E ’s need to be considered (this being 
impossible in practice, in e + e _ collisions the pro¬ 
cedure has an error of 0(a£ , ~ 1 ) if only ME’s with 
2 + n E < N partons are considered; in ref. [ 1551, 
N = 5). 


5.2. Complete matching of Monte Carlos 
and perturbative computations 

The problem of fully matching MC’s with 
higher-order computations can be seen as an up¬ 
grade of ME corrections: not only one wants to 
describe the kinematics of no + n E hard particles 
correctly, but the information on N raE LO rates 
must also be included. First attempts at solv¬ 
ing this problem have only recently become avail¬ 
able, and only for the case n E = 1. Following 
ref. [ 1611, let me denote by MC@NLO the im¬ 


proved MC we aim at constructing. The naive 
idea, of defining an MC@NLO by multiplying an 
MC with ME corrections by the NLO K-factor, 
is simply not acceptable: the inclusiveness of the 
K-factor does not really fit well into the exclu¬ 
sive framework of an MC. Thus, any MC@NLO 
must involve the computation of virtual ME’s. 
This is the reason why the construction of an 
MC@NLO is conceptually more complicated than 
ME corrections: the IR divergences of the vir¬ 
tual ME’s can only be cancelled by computing 
real-emission ME’s with one soft parton or two 
collinear partons. These soft and collinear con¬ 
figurations never occur in ordinary MC’s; in the 
case of ME corrections, the cuts S sep are specif¬ 
ically introduced to avoid them. The presence 
of virtual ME’s also requires a less-intuitive def¬ 


inition of double counting [ 161]: in the context 
of MC@NLO’s, double counting may correspond 
to either an excess or a deficit in the prediction. 
This generalization is necessary, since the ME’s 
used in MC@NLO’s are not positive-definite. 

In order to describe in more details current ap¬ 
proaches to MC@NLO, I adopt the toy model of 
ref. [ 161 , in which a system S (say, a quark) can 
emit “photons”, massless particles with only one 
degree of freedom (say, the energy). The initial 
energy of S' is 1, which becomes 1 — x after the 
emission of one photon of energy 0 < x < 1. The 
LO, virtual, and real contributions to the NLO 
cross section are (a is the coupling constant): 


da 

dx 


da\ 

dx A 


B5{ x), 

a(^+ v ) S(x), 


( 7 ) 

( 8 ) 
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da\ R(x) 

dx ) R x 


(9) 


where S(x) in eqs. (7 ; ) and (( 8 )) reminds that there’s 
no emission of real photons, and the IR divergence 
1 /e in eq. (||) results from the loop integration 
over virtual photon momentum in 4 — 2e dimen¬ 
sions. I denote by (5, z) the configuration of the 
system plus up to one photon, with z = 0 in the 
case of the LO or virtual contributions (eqs. (7i) 
and (||) respectively), and z = x ^ 0 in the 
case of the real contribution (eq. (j|)). Energy 
conservation is understood, and therefore in the 
configuration (5, z) the system has energy 1 — 2 . 
The function R(x) characterizes real emissions; 
its specific form is irrelevant, except that, for IR 
cancellation to occur, it must fulfil R{x) —> B for 
x —> 0. With one photon emission at most, any 
observable O can be represented by a function 
0(S , z)\ the computation of its expectation value 
(O) can be achieved through standard techniques 
for IR cancellation: 


(O) = £0(5,0) + . 


C B\og6 + V)0(S , 0) 


1-1 , 0(S,x)R(x) 
dx - 


( 10 ) 


in the slicing method [ , and 


(O) = f 

J o 


dx 


0(S, x) 


aR{x ) 


showers. In order to recover the correct total rate, 
each event resulting from a shower with initial 
condition (5, z) will be weighted with the coeffi¬ 
cient of 0(5, z) which appears in eq. ( 0 ) or in 
eq. ©• Using eq. (©, one gets 




Imc (O; 5, x) 


aR(x ) 
x 


I MC (0;S, 0) 


^B + aV — 



• ( 12 ) 


Unfortunately, this naive approach does not work. 
The weights aR(x)/x and B + aV — aB/x are 
IR-divergent at x = 0; since the corresponding 
showers have different initial conditions (S, x) and 
(5, 0), it would take an infinite amount of time to 
cancel the divergences (in other words, unweight¬ 
ing is impossible). This is not a practical prob¬ 
lem, is a fundamental one: the cancellation works 
for inclusive quantities, and the shower is exclu¬ 
sive. So the main problem in the construction 
of an MC@NLO can be reformulated as follows: 
how to achieve IR cancellation, without giving up 
the exclusive properties of the showers. Besides, 
eq. ( 0 ) also suffers from double counting. 

In the context of the slicing method, an ap¬ 


proach has been proposed [ |l62j, |163|, 164 , 
whic h exp loits an idea of refs. [ |166| 167 


165] 


see 


also ] |168| ). The slicing parameter 6 in eq. (10) is 
fixed to the value So, by imposing that no (5,0) 
contribution be present in the NLO cross section: 


+ 0(5,0) 


(b + clV 



( 11 ) 


in the subtraction method [ [kJ . 

In an MC approach, the system can undergo an 
arbitrary number of photon emissions. I denote 
by BI mc (0\S, 0) the distribution in the observ¬ 
able O obtained with MC methods; this notation 
reminds that in a standard MC the initial condi¬ 
tion for the shower is (5, 0) (the LO kinematics), 
and that the total rate is B (the LO rate, see 
eq. (7j)). The most straightforward implementa¬ 
tion of an MC@NLO can then be done by anal¬ 
ogy: since two kinematical configurations, (5, x) 
and (5,0), appear in the NLO cross section, one 
can use both of them as initial conditions for the 


B + a(BlogS 0 + V) = 0 . (13) 

This effectively restricts the energy of the real 
photons emitted to the range <5o < x < 1 
(see eq. ©)• This range is further partitioned 
by means of an arbitrary parameter (5 PS . One 
starts by generating the emission of a real photon 
with energy x distributed according to ciR(x)/x. 
Then, if Sq < x < S PS , the real-emission kine¬ 
matics (5, x) is mapped onto the LO kinematics 
(5, 0) (in other words, the photon with energy x 
is thrown away). The configuration (5, 0) is used 
as initial condition for the shower, requiring the 
shower to forbid photon energies larger than <5 PS . 
If x > S PS , the real emission is kept, and (5, x) 
is used as initial condition for the shower. The 
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corresponding formula is: 


% )=al ix 


' I mc (0-,S,x)^-0(x- 5 ps ) 


R(x) 


+/ MC (0; S, 0)-— 0(x - <5 O )0(<5 PS - x) 

x 


• ( 14 ) 


The advantage of eq. 0 ) is that it is manifestly 
positive-definite, and that its implementation can 
be carried out with little or no knowledge of the 
structure of the MC. On the other hand, it can 
be shown [ 161] that eq. 0 ) still has double 
counting, and does not have a formal perturba¬ 
tive expansion in a. These problems arise since 
the technique adopted to deal with the IR can¬ 
cellation is not exclusive enough: eq. (|l|) is an 
integral equation (the log <5o term is due to an in¬ 
tegral over soft-photon configurations in the real- 
emission contribution). Although the absence of 
a perturbative expansion can be seen as a minor 
drawback from a practical point of view, the im¬ 
pact of the double counting should be assessed 
for each observable studied, by considering the 
dependence of physic al p r edic t ions upon 5 PS . 

Other approaches [ 161, 16S| , 170, 171, 172, 173 
are based on the subtraction method. One can 
observe that ordinary MC’s do contain the infor¬ 
mation on the leading IR singular behaviour of 
NLO ME’s. Formally, this implies that the 0(a s ) 
term in the perturbative expansion of the MC re¬ 
sult can act as a local counterterm to the IR diver¬ 
gences at the NLO (strictly speaking, this is not 
exactly true in the case of large-angle soft gluon 
emission, and a few technical complications arise 
- see ref. [ 161 ). Furthermore, the form of the 
counterterm does not depend upon the observable 
studied. Thus, IR cancellation is achieved locally, 
but without any reference to a specific observable: 
this allows to implement it at the level of event 
generation, without giving up the exclusive treat¬ 
ment of the branchings. The prescription of ref. 


16l| is 


da 

dO 


= / dx 
Jo 


I MC (0;S, x) 


a[R(x) — BQ(x)] 


+7mc(0; S', 0) ( B + aV + 


aB[Q{x) - 1] 


where the Q(a;)-dependent quantities are the 
0(a) term of the MC result. The local IR can¬ 
cellation mentioned before shows in the fact that 
the coefficients of I M c(0; S, x) and I M c(0; S, 0) in 
eq. (|||) are finite, since the condition Q(x) —> 1 
for x —> 0 always holds, regardless of the spe¬ 
cific MC used. Therefore, these coefficients can 
be given as weights to the showers with (S, x) 
and (S, 0) initial conditions respectively. Eq. ( JTs] ) 
does not have double counting, and features a 
smooth matching between the soft- and hard- 
emission regions of the phase space, without the 
need to introduce any extra parameter such as 
6 PS . The price to pay for this is the presence 
of negative weights (which however do not spoil 
the probabilistic interpretation of the results). 
Furthermore, one needs to know details of the 
MC (Q(x)) in order to implement eq. ( 0 ): this 
seems unavoidable, since one should expect differ¬ 
ent MC’s to match differently with a given NLO 
computation. The first QCD implementation of 
eq. 0) has been presented in ref. 


The approach of refs. [ |l69[ |17C 


a technique similar to that of ref. 
on the definition of a local IR counterterm. At 


] 174 

I- 


171, 

172 

uses 

] 16 

1 ], based 


variance with ref. [ 161], where the resummation 


of large logs is performed to LL accuracy, refs. 
169, 170, 171, 172] advocate an NLL (or beyond) 
resummation; for this to happen, it is argued that 
the standard formulation of collinear factoriza¬ 
tion must be extended. This approach has been 


fully formulated [ 1711 only in the unphysical </>]) =6 
theory so far. Current QCD implementations do 
not include gluon emission. 

6. CONCLUDING REMARKS 

It seems appropriate to start this section by 
mentioning a couple of phenomenological is¬ 
sues which would have deserved more attention. 
One is the problem affecting the single-inclusive, 
isolated-photon measurements at the Tevatron. If 
DO data are considered f 175, 176], a moderate 


(15) 


disagreement with NLO QCD is present in the 
low-p T , central- r] region, which disappears when 
the ratio R = a[\fS = 1800 GeV) / a(\fS = 
630 GeV) is considered. The situation worsens 
in the case of CDF data [ 177]: not only the dis- 






















































16 


crepancy with QCD is statistically more signif¬ 
icant for cross sections, but the measured ratio 
R also disagrees with theory. The consistent in¬ 
clusion of recoil effects, along the lines of refs. 
[ 178, 176], might increase the QCD prediction 
at small p T ’ s, and thus reduce the discrepancy. 
The second problem affects the single-inclusive jet 
cross section as measured by DO [ 18C], the jets 
being reconstructed with a fc T -algorithm. The 
data display a rather poor agreement with the¬ 
ory for p T < 100 GeV; this is disappointing, given 
the excellent results obtained with the cone algo¬ 
rithm (as far as the cone algorithm is concerned, 
I should also mention here that the previously 
reported excess of data over theory at large p T 
has now completely disappeared: NLO QCD per¬ 
fectly reproduces the data, if updated PDF sets 
are used. An explanation in terms of PDFs has 
been already given in the past, but the PDF set 
used at that time resulted from an ad hoc fit - 
named “HJ” by the CTEQ collaboration. This is 
not necessary any longer, since in the newest PDF 
releases the gluon density of the best fit naturally 
results to be HJ-like. See ref. [ 181 for a dis¬ 
cussion on this point). It is hard to believe that 
the discrepancy in the case of the fc T -algorithm 
is the signal of a serious problem in QCD (since 
it would probably affect the cone algorithm as 
well); however, it may indicate a deficiency in 
current MC simulations, or in the understanding 
of lradronization and/or detector effects, which 
would surely worsen at the LHC energies. It has 
to be remarked that the fc T -algorithm appears to 
work well at HERA. 

In general, the capability of QCD to de¬ 
scribe hard production processes is quite remark¬ 
able (with the exception of 6-production at LEP 
and HERA). The definition of a formalism for 
the computation of exclusive observables at the 
NNLO is one of the most challenging and hot 
topics in perturbative QCD, which will have an 
important impact on phenomenological studies 
and will be crucial in improving the precision 
of a s (and other fundamental parameters) mea¬ 
surements - at colliders, exclusive Drell-Yan pro¬ 
duction will necessarily have to be calculated to 
NNLO, in order to match the experimental pre¬ 
cision. The substantial progress made in the past 


couple of years in the computation of two-loop in¬ 
tegrals and three-loop splitting functions, and the 
NNLO results for direct inclusive SM-Higgs pro¬ 
duction, are certainly very encouraging. NNLO 
results will also help to understand better the in¬ 
terplay between soft and hard physics; the e + e~ 
and DIS environments will serve as a laboratory 
for the more involved case of hadronic collisions, 
where most of the work remains to be done. 
Monte Carlos are much better equipped to face 
the challenges of Tevatron Run II and LHC. The 
implementation of new processes in matrix ele¬ 
ment generators proceeds steadily and, although 
there is still work to be done on the general struc¬ 
ture of matrix element corrections, these tech¬ 
niques are rather well established by now. For¬ 
malisms for MC@NLO’s have received consider¬ 
able attention in the recent past; the field has still 
to reach a mature stage, and new ideas will cer¬ 
tainly be presented in the near future. A rather 
obvious but quite challenging step is that of merg¬ 
ing matrix element corrections and MC@NLO ap¬ 
proaches. 

It is a pleasure to thank my friends and col¬ 
leagues, unfortunately too numerous to mention, 
who helped me so much in the preparation of the 
talk and of this manuscript. 
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